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Abstract. We present cosmological simulations of galaxy clusters, with focus on the clus- 
ter outskirts. We show that large-scale cosmic accretion and mergers produce significant 
internal gas motions and inhomogeneous gas distribution ("dumpiness") in the intracluster 
medium (ICM) and introduce biases in measurements of the ICM profiles and the cluster 
mass. We also show that non-thermal pressure provided by the gas motions is one of the 
dominant sources of theoretical uncertainties in cosmic microwave background secondary 
anisotropies. We briefly discuss implications for cluster cosmology and future prospects 
for understanding the physics of cluster outskirts using computer simulations and multi- 
wavelength cluster surveys. 



1. Introduction 

In recent years, galaxy clusters have emerged 
as one of the most unique and power- 
ful laboratories for cosmology and astro- 
physics. Being the largest and most mag- 
nificent structures in the Universe, clusters 
of galaxies serve as excellent tracers of the 
growth of cosmic structures. Current genera- 
tion of X-ray cluster surveys have provided 
independent confirmation of cosmic acceler- 
ation and significantly tig hten constra i nts on 
the nature of dark en ergy dAllen et all l2008t 
Vikhl inin et all l2009l) and alternative theories 
of gravity (e.g., ISchmidt et all l2009h . Several 
ongoing and new X-ray (e.g., eROSITA) and 
Sunyaev-Zel'dovich effect (SZE) cluster sur- 
veys (e.g., SPT, ACT, Planck) are underway to 
improve current cosmological constraints. 

Outskirts of galaxy clusters have special 
importance for cluster cosmology, because 
they are believed to be much less suscepti- 
ble to complicated cluster astrophysics, such 



as radiative gas cooling, star formation, and 
energy injection from active galactic nuclei. 
Dominant physical processes in the outskirts 
are limited to the gravity-driven collisionless 
dynamics of dark matter and hydrodynamics 
of the intracluster medium (ICM). In the hi- 
erarchical structure formation model, galaxy 
clusters grow by accreting clumps and diffuse 
gas from the surrounding large-scale struc- 
ture in their outer envelope. Numerical simu- 
lations predict that the large-scale cosmic ac- 
cretion and mergers give rise to internal gas 
motions and inhomogeneous gas distribution 
in the ICM. However, until very recently, ob- 
servational studies of the ICM have been lim- 
ited to radii considerably smaller than the virial 
radius of clusters. 

Recently, Suzaku X-ray observations have 
extended X-ray measurements of the ICM pro- 
file out to an d beyond the virial radius for sev- 
eral c l usters dBautz et al. , 2009); iGeorge et ail 
2009; Reim-ic h et all l2009t Hoshino et all 
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l2010t iKawaharada et all 1201 Oh . While these 
measurements are still quite uncertain, initial 
results suggested that the observed ICM pro- 
files may deviate significantly from the pre- 
dictio n of hydrodynamica l cluster simulations 
(e.g., iGeorge et all [2009). In addition to test- 
ing models of structure formation, these new 
measurements will be important for controlling 
systematic uncertainties in cluster-based cos- 
mological measurements. 

In this work, I will present theoretical 
modeling of the outskirts of galaxy clusters 
based on cosmological simulations, with high- 
lights on implications for the interpretation of 
forthcoming multi-wavelength observations of 
galaxy clusters. The simulations we pre sent 
here are described in Naga i et alJ (|2007a) and 
INagai et ail (l2007bl) . and we refer the readers 
to these papers for more details. 

2. Gas Motions in Cluster Outskirts 

Gas motions induced by cosmic accretion and 
mergers provide non-thermal pressure support 
in galaxy clusters. The top panel of Fig. [2] 
shows that results of hydrodynamical cluster 
simulations, illustrating that the non-thermal 
pressure contributes to 10-20% at r — rsoo of 
the total pressure. The fraction of non-thermal 
pressure support increases with radius, and it 
is larger for more dynamically active systems. 
The bottom panel of Fig. [2] shows that the 
cluster mass profile based on hydrostatic as- 
sumption is biased low. We further demon- 
strate that it is possible correct the bias in the 
hydrostatic mass, if one could measure gas 
motions in clusters and hence correct for the 
bias. Upcoming Astro-H X-ray satellite mis- 
sion (scheduled to be launched in 2014) will 
provide a first direct measurement of gas mo- 
tions i n clusters via doppler broade ning of iron 
lines (llnogamov & SunvaevLl2003l) . 

3. Gas Clumping in Cluster Outskirts 

Recently, Suzaku X-ray observations revealed 
that the observed entropy profile of the ICM is 
significantly offset from the prediction of hy- 
drodynamical simulations of galaxy clusters. 
Here, we point out that gas clumping is likely 
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Fig. 1. Top panel: Ratio of pressure from random 
gas motions to total pressure as a function of radius. 
Relaxed clusters are represented by solid lines while 
unrelaxed clusters are represented by dashed lines. 
Bottom panel: Averaged mass profiles M{< r) of the 
relaxed clusters, normalized by Msoa- The solid line 
shows the actual mass profile from simulation, the 
long dashed line shows the mass profile from hydro- 
static equilibrium including random gas and thermal 
pressure, and the short dashed line shows the mass 
profile from hydrostatic equilibrium including ther- 
mal pressure only. Ha shed region show s the 1-cr er- 
ror of the mean. From lLau et ail ([2009). 

a major source of systematic bias in X-ray 
measurements of ICM profiles i n the envelope 
of ga laxy clusters (r > r2oo) (INagai & LauL 
2011). Using hydrodynamical simulations of 
cluster formation, we show that gas clumping 
introduces the overestimate of the observed gas 
density and causes flattening of the entropy 
profile at large radius. This is illustrated in 
Fig. [2] The top panel shows that the clumping 
factor of the X-ray emitting gas (T > 10 6 K) 
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Fig. 2. Top panel: Median clumping factor pro- 
files of gas with different minimum temperature for 
the simulations with gas cooling and star forma- 
tion. Bottom panel: Impact of gas clumping on X- 
ray measurements of the ICM entropy profile. The 
dashed line indicates the true entropy profile of the 
simulated clusters, while the solid line indicates the 
observed entropy profile, obtained by assuming no 
clumping (C = 1). Black points are Suzaku observa- 
tions of PKS0745-191 (circles), A1689 (triangles), 
A141 3 (stars), and A1795 (black dot-dashed lines). 
From lNagai & Laul d201ll) . 



is C = <Pg as )/(pg as ) 2 ~ 1.3 at r = r 20 o, 
and it increases with radius, reaching C ~ 5 
at r — 2r2oo- In the bottom panel, the solid 
line indicates the true entropy profile, which 
is consistent with the self-similar prediction, 
K = T/n 2 e ' 3 cc r u dVoit et all 120051) . From 
the definition of entropy, the overestimate of 
gas density due to clumping causes an under- 
estimate of the observed entropy profile by 
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Fig. 3. Astrophysical uncertainties in the SZ power 
spectrum. Different lines indicate theoretical uncer- 
tainties associated with (1) heating of gas by energy 
injection from stars and AGN (indicated with solid 
lines) and merger dynamics (indicated with dashed 
lines in top-left panel), (2) dark matter structures 
(top-right), and (3) non-thermal pressure by gas mo- 
tions (bottom-left) and its time evolution (bottom- 
right). In each case, t he thick red line re presents our 
fiducial model. From Shaw et al. 



C(r) 1 / 3 . Results of our analyses indicate that 
gas clumping causes the flattening of the ob- 
served entropy profiles at r > r^. While cur- 
rent Suzaku measurements are still uncertain, 
our results indicate that gas clumping is im- 
portant for reducing the tension between recent 
Suzaku observations and theoretical prediction 
of the ACDM model. 

4. Impact on the SZ power spectrum 

Recent measurements of the SZ power spec- 
trum by SPT and ACT telescopes revealed 
that the SZ power is significantly below the 
signal predicted by the c urrent cosmi c struc- 
ture formation model dLueker et all 120101: 
IShirokoff et aill201QtlDunklev et al.Ll2010h . 

In our recent work, we argued that the cur- 
rent SZ power spectrum template is overesti- 
mated by 50-100%, due to lack of important 
astrophysical process in th eoretical modelin g 
of the SZ power spectrum (IShaw et al.L l2010h . 
Fig. [3] illustrates theoretical uncertainties in 
the thermal SZ power spectrum. Our model 
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accounts for star formation and energy feed- 
back (from supernovae and active galactic nu- 
clei) as well as radially dependent non-thermal 
pressure support due to random gas motions, 
the latter calibrated by recent hydrodynami- 
cal simulations. Varying the levels of feedback 
and non-thermal pressure support can signif- 
icantly change both the amplitude and shape 
of the thermal SZ power spectrum. Increasing 
the feedback suppresses power at small angu- 
lar scales, shifting the peak of the power spec- 
trum to lower I. On the other hand, increas- 
ing the non-thermal pressure support signifi- 
cantly reduces power at large angular scales. 
In general, including non-thermal pressure at 
the level measured in simulations has a large 
effect on the power spectrum, reducing the am- 
plitude by > 60% at angular scales of a few ar- 
cminutes compared to a model without a non- 
thermal component. Comparing with the re- 
cent measurements of the small-scale cosmic 
microwave background power spectrum, our 
model reduces the tension between the values 
of <t 8 measured from the SZ power spectrum 
and from cluster abundances. 

5. Future Prospects 

Modern numerical simulations predict that gas 
motions and clumping are ubiquitous in the 
outskirts of galaxy clusters. While they are 
generic predictions of the concordance ACDM 
model, we have had very little observational 
handle on these phenomenon until very re- 
cently. 

New X-ray and SZE observations just com- 
ing online have significantly extended mea- 
surements of the ICM, out to and beyond the 
virial radius of clusters. For example, deep X- 
ray imaging of nearby clusters with current 
Suzaku and Swift-XKY have started to pro- 
vide accurate measurements of the ICM pro- 
files in the outer envelope of galaxy clusters. 
Comparison of current X-ray and SZE mea- 
surements should soon provide insights into 
the properties of the ICM in cluster outskirts. 
Astro-H X-ray mission promises to provide the 
first direct measurements of internal gas mo- 
tions in clusters, and eROSITA and WXRT X- 
ray mission will help increase the number of 



clusters with deep imaging data extending out 
to large radii. Theoretical modeling based on 
detailed numerical simulations is also under- 
way. A plethora of activities (in both theory 
and observation) will advance our understand- 
ing of cluster physics and provide foundation 
for the use of galaxy clusters as laboratories for 
cosmology and astrophysics. 
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